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Abstract 

The geometric features of the square and triadic Koch snowflake drums are compared 
using a position entropy defined on the grid points of the discretizations (pre-fractals) 
of the two domains. Weighted graphs using the geometric quantities are created and 
random walks on the two pre-fractals are performed. The aim is to understand if the 
existence of narrow channels in the domain may cause the 'localization' of eigenfunc- 
tions. 



1 Introduction 



The term fractal is defined very loosely [Fj. Here are just some ways of characterizing a 
fractal set F: 

• F has a fine structure, i.e., detail on an arbitrary small scale. 

• F is too irregular to be described in traditional geometric language, both locally and 
globally. 

• Often, F has some form of self-similarity (maybe strict, approximate or statistical). 

• Usually, the fractal dimension of F (defined in some way) is greater than its topological 
dimension. 

• In most cases, F is defined in a very simple way (e.g., recursively). 

In short, fractals are objects with irregular geometry. Such objects can be found every- 
where in nature, the most obvious example being that of a tree structure. The vibrational 
properties of fractals are of great interest. Questions that arise naturally are, for example: 

• Why are waves damped much more by fractal rather than by smooth coastlines? 

• How can we explain the vibrational properties of glass? 

• Why does fractal foliage of trees provide such strong resistance to the wind? 

All these questions are largely unanswered. Fractal geometry is a means to describing 
strongly irregular objects and questions like the ones mentioned above have played and 
continue to play an important role in the development of the subject. The main idea is to 
use strongly irregular but deterministic objects as a good approximation to the fractal. If 
the physical properties of the object of study are related to the character of the geometry, 
these deterministic approximations prove to be a good source of information. 



1.1 Fractal Drums 

A fractal drum is the simplest example of a surface fractal resonator. It is a flat surface 
bounded by a fractal frontier. 



Figure 1: An approximation to the triadic Koch Snowflake domain 



1.2 Vibrational properties of fractal drums 

The vibrational properties of these structures are of great theoretical and practical interest 
HEH, [Big, |Hro(^a| . (Ll], |L2|, jOl, (LHESI, |LP|, [E^RPSJ. |SII1, [^HM] . 

An extensive amount of experimental work on the study of waves or harmonic oscillations 
carried by fractals has been done by the physicist B. Sapoval and his collaborators |SG| . 
|S(tM| . The fractal drum in their experiments consisted of a stainless steel sheet with the 
boundary etched in the pattern of a square Snowflake. A soap bubble deposited on this 
object was excited to low frequency vibrations by a loudspeaker situated above the fractal 
drum. Their results showed that the wave motion is strongly damped by the fractal bound- 
ary and that localization may occur. Sapoval attributed this localization to the existence of 
narrow channels in the geometry of the square domain. This phenomenon was not observed 
on the triadic Snowflake domain. 

In Snowflake harmonics and computer graphics: numerical computation of spectra on fractal 
drums |LNR,G| . the authors numerically estimated the first 50 eigenvalues and eigenvectors 
of the triadic Snowfiake domain and tested graphically the results and conjectures concerning 
the boundary behavior of the gradient of the eigenfunctions. 

Some mathematical work has been done to confirm Sapoval's results. Recently, M. L. Lapidus 
and M. Pang |LPj have studied the boundary behavior of the Dirichlet Laplacian eigenfunc- 
tions and their gradients on a class of planar domains with fractal boundary, including the 
(triangular and square) Koch Snowfiake domains, as well as their polygonal approximations. 
One of their main results, specialized to the Koch Snowfiake domain, states that the mag- 
nitude of the gradient of the first eigenfunction (or 'first harmonic') 'blows up' at infinitely 
many boundary points, i.e., the membrane of the Koch snowfiake drum exhibits 'infinite 
stress' near such points. Physically, this corresponds to a strong damping phenomenon. 
In this paper, we develop a new approach to investigate localization phenomena on fractal 
drums: the discretizations of the square and triadic Snowfiake are considered as weighted 
graphs. The weight is assigned to an edge in such a way that it refiects certain geometric 



properties of the two nodes connected by that edge. The development of random walks 
initiated in different geometric regions are compared and analyzed. 



2 Measures of localization 



The main problem when studying the localization of the eigenmodes is to define a suitable 
notion of localization. What do we mean when talking about localization of an eigenfunction? 
A function on an unbounded domain is said to be localized if it has compact support or if it is 
exponentially localized but a generally accepted notion of localization on bounded domains 
is still missing. In fact, we note that due to the ellipticity of the underlying equation, there 
cannot exist any Dirichlet eigenfunction that vanishes on a nonempty open subset of the 
domain |Evj . 

2.1 The square and triangular Koch prefractals 

The square and triangular prefractals are defined recursively by repeatedly applying the cor- 
responding generators. It was proved that the sequence of square and triangular prefractals 
converge to the corresponding fractals and that the sequences of the individual eigenvalues 
of the prefractals converge to the eigenvalues of the actual fractals. Thus, the study of the 
eigenfunctions on these approximations should give a good insight into the behavior of the 
eigenfunctions on the true fractals. 

Figure 2: Square Snowfiake Prefractal 



Figure 3: Triadic Snowflake Prefractal 




2.2 Localization measures 

Some eigenmodes on the square Koch prefractal show confinement to small regions of the 
membrane. On triadic Koch prefractals, however, no such 'localization' is observed. Sapoval 
attributed the experimental localization he observed to the effects of 

1. Damping 

2. The existence of 'narrow paths' in the geometry of the structure. 



Figure 4: Eigenfunction 27 



The following quantities were introduced in the literature in an attempt to qualitatively 
describe localization of functions (^) on bounded domains: 



1 ) Participation Ratio Vp(\E') [WeJ 
The Participation Ratio is defined as PR{^) 



In I'fi^Tdx 



1 



. It can be regarded as the effective 



number of nodes participating in the function \i/ with significant weight. 

2) Participation Volume ^(^') |BeDe] 
The Participation Volume is a measure of the effective number of 'atoms' of an 'ensemble' 



the volume of the domain D, can be used to distinguish between 'extended' (non-localized 
states) and 'exponentially localized states' (the amplitude decays exponentially for suf- 
ficiently large distances from some central point). 

In this paper we use the Participation Ratio in conjunction with suitable notions of diame- 
ter and entropy to investigate localization phenomena on the square and triadic Snowflake 
domains. 

3 The Geometry of the Square and Triadic Snowflake 
Domains 

How can we describe the existence of narrow channels in the geometry of the square Snowflake? 
We will attempt to classify different regions by defining the following six distances for an 
interior grid point {x,y): 

Let dh^{x, y), dh_{x, y) be the right and left horizontal distances of (x, y) from the boundary 
and (it>+, dv_ be the up and down vertical distances of (x, y) from the boundary. Now define 

• dh{x,y) = dh+{x,y) + dh^{x,y) 

• dv{x,y) = dv+{x,y) + dv_{x,y). 



participating in a vibration on a domain D: Vp(\E') 



Vp(^)/V, where V is 



Using dh{x, y) and dv{x, y), one can distinguish between three geometrically different regions 
of the domain in question: 

1. 'The Grottos': Regions confined in both directions. 

For example, some gridpoints close to the boundary have dv ~ dh taking small values. 

2. 'The Canyons': Narrow channels. 

Points lying in those regions have either dv <^ dh or dh <^ dv. 

3. 'The Prairies': Relatively unconfined regions. 

Some points in the center of the domain, for example, have dv,dh ~ 1/2 and hence 
should be considered as lying in an unconfined region. 

Remark: Note that points in region one as well as points in region three have dv ~ dh. 
Since the ultimate goal is to define a measure distinguishing between the three different 
geometric regions, it is clear that plainly using dh{x,y) and dv{x,y) will not do. 

We thus define the following three quantities: 

1. dmin{x, y) = min {dh+{x, y), dh_{x, y), dv+{x, y), dv_{x, y)} 
z. Urat{3[,, y) — mm \ ^^^^ ^^ , ^^^^^^^ j- 

3. (x, y) = drat{x, y) * dmin{x, y). 

Note that for interior grid points in SQ{L,R), (L > 0) with grid size Hsq(l,r) — (|)^^^^: 
we have 

1- djnin e [HsQ{L,R)A/^] 

2. drat e [2HsQ(L,R), 1] 

3. a-^e[2(i^sQ(L,«))^l/4]. 



For interior grid points in T{L, R), {L > 0) with grid size Ht{l,r) — {\) , we have 

2. drat e [-^Ht(L,R), 1] 

3. a-^e[(i^T(L,i^))^^]. 

Let us consider the properties of drat and d^!^^ in the three different regions: 
Properties of drat 

• if {x, y) lies in a grotto, 

then dv dh and drat{x, y) ~ 1. 

• if {x, y) lies in a canyon, 

then dv <^ dh or dh <^ dv and drat{x, y) is small. 

• if {x, y) lies in a prairie, 

then dv dh and drat{x, y) ~ 1. 

Properties of 9^*" 

Let us assume that the fractal and refinement levels are not both zero. 

• if {x, y) lies in a grotto, 

then dmin{x,y)) is small. 

• if (x, y) lies in a canyon, 

then dl!X{x,y) (f^ {dmin{x,y)f) is very small. 

• if {x, y) lies in a prairie, 

then dli!^^{x,y) (~ dmin{x,y)) is relatively large. 



Let us now compare the features of y) on the two domains: 



Table 1: Comparison of a™"(a;,|/) at {L,R) = (2,2) 




Square Snowflake 


Triadic Snowflake 


Maximum 


0.25 


0.4208 


Minimum 


0.001 


0.0002 


Minimum/Maximum 


4 * 10-=^ 


4.75 * 10"^ 


Median 


0.0156 


0.0216 


Mean 


0.0332 


0.0724 


Standard Deviation 


0.0467 


0.0967 


Median /Maximum 


0.0624 


0.0513 


Mean/ Maximum 


0.1328 


0.172 



Note from table 1 that the maximum of d^!^^ on the square Snowflake is | while it is 
on the triadic Snowflake. The minimum value of 9™^" on the interior of the triadic domains 
is smaller than that on the square domain, while the mean value of c?™j" is smaller on the 
square Snowflake. 

We find that in both domains, the mean value of 9^^" lies closer to the minimum value 
of 9^^" than to the maximum value. Together with the low medians (only 6.24% of the 
maximum on the square, 5.13% of the maximum on the triadic domain), this means that a 
lot of points on either domain have a relatively small 9™/' value. The relative mean value 
of d^!^^ is higher on the triadic Snowflake due to the triangular geometry. At reflnement 
and fractal level equal to two, we flnd that the difference between the maximum and the 
mean value lies at about 0.2168 for the square Snowflake and at about 0.4206 for the triadic 
Snowflake. The differences between minimum value and mean value lie at about 0.0322 and 
0.0722, respectively. 



Figure 5: Square Snowflake c?™" distribution 
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Figure 6: Square Snowflake d^^J^ surface plot 



Figure 7: Triadic Snowflake dJ^^^ distribution 




Figure 8: Triadic Snowflake dJ^^^ surface plot 




4 Position Entropy 

In this section, we will use 5™^", defined in the previous section, to define a probability 
distribution P^^yid) on each node of our discretization. This probability distribution is used 
to define a position entropy S{x, y) at each vertex. 



For {x,y), a grid point of the discretization of the snowflake domain, let 

P(x,y){9) '■— percentage of neighbours of node {x,y) with (9™™ = 9. (1) 

Then '^gP{x,y){d) — 1 for each node {x,y) and hence, P defines a probability distribution. 
Now let 

S{x, y):^-J2 Pi^,y)(9) In (2) 
a 

This entropy gives information about the disorder in the values of 9^^" at neighboring points 
of the vertex in question. The more variation in 9^^", the larger the value of S{x, y) will be. 



Table 2: Comparison of Position Entropy at (L, R) = (2, 2) 



Position Entropy 


Square Snowflake 


Triadic Snowflake 


Maximum 


0.0157 


0.0238 


Minimum 


2.63 * 10-1^ 


1.28* 10-1° 


Minimum/Maximum 


1.675 * 10-12 


5.378* 10-9 


Median 


0.0085 


0.0114 


Mean 


0.0089 


0.0114 


Standard Deviation 


0.0041 


0.0059 


Median/Maximum 


0.5414 


0.47899 


Mean/Maximum 


0.5668 


0.47899 



Comparing our two domains at {L,R) — (2,2), we find from table 2 that tlie mean and 
median values are the same on the triadic Snowflake and close to each other on the square 
Snowflake, and in both domains, they lie closer to the maximum value of S than to the 
minimum value. This means that most of the points on either domain have a relatively large 
entropy values, i.e., the values of dj^^ in the neighborhood of those points are not evenly 
distributed. The relative median and mean of S are higher on the square Snowflake. At 
refinement and fractal level equal to 2, we find that the difference between the maximum 
and the mean value lies at about 0.28 for the square Snowfiake and at about 0.38 for the 
triadic Snowfiake. The differences between minimum value and mean value lie at about 0.83 
and 0.95, respectively. 

In summary, we can say that the relative mean and median of the entropy found on the square 
domain are significantly larger ( about 10%). The d^!^^ characteristics of both domains are 
similar but we find a slightly higher (1%) relative median and a lower mean ( about 4%) on 
the square Snowfiake. 

We refer to figures 9-12 for the distribution and the surface plot of the square and triadic 
Snowfiake domains. 



Figure 9: Square Snowflake Position Entropy distribution 
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Figure 10: Square Snowflake Position Entropy surface plot 



Figure 11: Triadic Snowflake Position Entropy distribution 

140 I 1 1 1 1 1 



120 - 



100 - 



c 
o 




0.005 0.01 0.015 0.02 0.025 

Position Entropy 



Figure 12: Triadic Snowflake Position Entropy surface plot 



5 Random Walks on Discretizations of Snowflake Do- 
mains 



In this section, we will study diffusion (random walks) on discretization graphs of the 
Snowflake domains. This work is inspired by I. Simonsen et. al. |SiEr MaSn] . [Si) . In 
their papers, these authors use diffusion on networks to derive a master equation whose 
analysis reveals information about the large-scale structure of the network. By performing 
random walks governed by the underlying master equation on the two Snowflake domains, 
we hope to gain useful insights into their geometric properties and an understanding of why 
localization is favored by the Square domain. This will provide a new way to investigate 
localization for domains with fractal boundaries. This study relies in part on the notions of 
distance and entropy introduced earlier in this paper (see Sections 3 and 4). 

The nodes of the discretizations represent the vertices and the edges are given by the neighbor 
relations. The path taken by our random walkers will be influenced by the geometric features 
of the domain: 

The edge Cij, going from vertex Vi to vertex vj is assigned a weight inversely proportional to 
1 — <9^™(j), where (9^™ is normalized such that it takes values in [0, 1] . Hence, edges where 
node Vj lies in a Grotto have large weights assigned, edges where node Vj lies in a Canyon 
have even larger weights (close to 1) assigned, while edges with node Vj lying in a Prairie have 
relatively small weights assigned, i.e., the weights are assigned to the edges in such a way 
that a random walker is strongly encouraged to move to regions with distances that describe 
narrow channels, then to regions confined in both directions, and lastly to unconfined regions. 
If localization is caused by the existence of narrow channels, as predicted by Sapoval, we 
expect to see some kind of congregation of the walkers in those areas of the domain that are 
characterized by predominantly large weights. 

We start our random walk by placing a large number of random walkers onto the vertices 
of the network. At each time step, these walkers are allowed to move between adjacent 



vertices. The edge, out of the possible outgoing ones, a walker chooses to move along is 
picked at random with a probability equal to the weight assigned to this edge. 



We define the adjacency matrix A as follows: 



1, if nodes Vi and Vi are neighbors 
0, else. 
Note that the weight matrix W — {wij}ij = — ^rar(i))} non-symmetric. 

5.1 The Master Equation 

Let nv be the number of vertices in the discretization and be the number of participating 
random walkers. Let the number of walkers situated at node Vi at time t be denoted by 
Ni{t). Then the fraction of walkers at this node at time t is 

,m = (3) 
Since the total number of random walkers is preserved at each time, we have ^j?7i(t) = 1. 

The change in the walker density of a vertex Vi during one time step equals the difference 
between the relative number of walkers entering and leaving the same vertex over the time 
interval. In mathematical terms we can write 



rj,{t+l)^rjiit) + Jrit)-J^{t), (4) 
where J^{t) denote the relative number of walkers entering (-) and leaving (+) vertex Vi. 
The value of Jj^{t) depends on the strength Sti of node f » i.e., the total number of outgoing 

weights Sti — ^^j^Neighj 

Wij, where Neighi denotes the set of neighbors of node Vi. Similarly, 
the value of Jf{t) depends on the total number of incoming weights YlijeNeighi'^ji- '^^^ 



fraction of outgoing walkers from vertex Vi (a current) per unit weight is thus 

= . (5) 

Hence, Cj(t) is the weighted walker density per link of node Vi at time t. 

The edge current on the directed edge from vertex Vi to vertex Vj is then given by 

Cijit) = WijCiit) = Wij— . (6) 



Notice that the term 



is the probability of a walker moving from vertex Vi to vertex Vj . 

The relative number of outgoing walkers from node Vi at time t is given by 



and the relative number J~ of incoming walkers to node Vi at time t is given by 



J-(t)= (8) 

j^Neighi 

It is easy to verify that J^(t) = rjiit). This expresses the fact that all walkers leave their 
respective vertices at each time step, i.e., no walker stays at the same vertex for more than 
one time step. 

Denoting dtr]i{t) = rjiit — 1) — r]i{t) and substituting the expressions for J^{t) and J~(t) into 
equation (jH), we find 



j&Neighi 



(9) 



j^Neighi 



where 



l^kf^Neighj ^jk 

This equation can be written in matrix form, as follows: 



(10) 



Equation (fTfl|l is known as the master equation for the random walk process of the under- 
lying network. Further, T = {Tij} is called the transfer matrix. It 'transfers' the walker 
distribution one step ahead and therefore can be thought of as a time propagator for the pro- 
cess. For some arbitrarily chosen initial state ^^(O), the time development can be obtained by 
iterations of equation JTH)), with the result r]{t) = T^^^r](0). Hence, the eigenvalue spectrum 
of T controls the time evolution of the diffusive process. 

Due to the fact that our weight matrix is not symmetric, the matrix T also fails to be 
symmetric. However, T is similar to the symmetric matrix S = KTK^"*^, where 




and hence 




y/{i-dzrmi-dzr{j)) 



(11) 



Therefore, T is guaranteed to have real eigenvalues and corresponding eigenvectors. It 
is convenient to sort the eigenvalues in descending order. Furthermore, as a consequence 
of the Perron-Frobenius theory, the eigenvalue Ai = 1 is simple and the elements of the 
corresponding eigenvector will all have the same sign. 

Physically, the state Ai = 1 corresponds to the stationary state such that 77(00) oc 77*, where 
the diffusive current flowing from node Vi to node Vj is exactly balanced by that flowing 
from Vj to Vi. Since |Afe| < 1 for A; 7^ 1, all modes corresponding to those eigenvalues are 
decaying. A^ > correspond to non-oscillating modes while A^ < correspond to states 
where oscillation will take place over time. The large scale topology of the given network 
reflects itself in the statistical properties of the eigenvectors 77^ fS iErMaSn| . 

6 Numerical Work 

As shown in the previous section, a random walk on the discretization of the square and 
triadic Snowflake domains is governed by the master equation ri{t) = T^^^t](0), where r](t) is 
the distribution of random walkers at time t. In this section, we will discuss numerical sim- 
ulations of these random walks. Initially, we will place all random walkers at one vertex. By 
repeatedly applying T to our distribution, we are able to obtain the new walker distributions 
for each time step. We expect that the dynamics of this system is strongly influenced by 
the choice of the initial node. Recall that a node lying in a narrow channel is characterized 
by small 9^^*^ and entropy values. Hence, in accordance with the definition of the weight 
matrix W {wij = 1 — drat^in{j)) , if narrow channels in the geometry of the square Snowflake 
exist, we would hope to see a concentration of walkers in the regions where 9^^"^ and entropy 
values are low. 



Figure 13: Triadic Snowflake: Position Entropy plus Ratio Distance contour plot 




In figure 13, we see a plot of the relative Minimum Ratio Distances added to the relative 
Position Entropies at each grid point of the Triadic Snowflake. As expected, we see low 
values close to the boundary. Indeed, random walkers situated here are expected to stay 
localized in the region. The center part of the Snowflake comprises a medium value square 
on the left and right of which we find two rectangular discrete low bands and 4 other medium 
value regions situated north, west, south and east of the center square. In summary, the 
triadic Snowfiake has a relatively big open area of low and medium values in its center. This 
is not the case for the square Snowfiake. 



Figure 14: Square Snowflake: Position Entropy plus Ratio Distance contour plot 




Figure 14 shows a plot of the relative Minimum Ratio Distances added to the relative Position 
Entropies at each grid point of the square Snowflake. Note that in the center we find relatively 
high values but, surrounding that center, we find a band of low entropy and distance values. 
We see four arms of high entropy/distance stretching through parts of the domain almost 
reaching the boundary. In each of those four arms we can find a small low value basin. 
The remaining regions between the arms and near the boundary are also characterized by 
low values. We would expect that a random walk started in any of the basins would stay 
localized, i.e., that the majority of high/low amplitudes can be found there. If the random 
walk is initiated in a high altitude region, we would hope to see no such localization. If this 
were observed, it would indeed indicate that localization really is caused by the existence of 
narrow channels. 



To test this theory, we initially place all random walkers at a point of interest (low en- 
tropy/low distance, high entropy/high distance, low entropy/high distance or high entropy/low 
distance) and then watch the random walk evolve by computing the following quantities at 
time t: 

• the diameter DIAM(t), the furthest distance from the initial node traveled by a ran- 
dom walker during [0,t]. 

• the Participation Ratio PR{t) = ^ [^ (t)]'^ ' where r]i{t) is the percentage of random 
walkers at node Vi at time t. 

We define high amplitude points at time t as points where \r]i{t)\ > {3/4:)'max{\r]i{t)\) and 
compute 

• D'^j^^it), the Relative Minimum Ratio Distance of high amplitude points at time t, i.e., 
Dhtit) = "i^f jgHn)"^*^^ ^ where mean((9;^f is the mean of the Minimum Ratio 
Distances of high amplitude points. 



• ^hai^)^ ^he Relative Position Entropy of high amplitude points at time t, i.e., S'f^^{t) — 
'^^°max(s)^^ , where mean{Sha{t)) is the mean of the entropies of high amplitude points. 

• DIAM^l^t) — ^^j^j^i the Relative Diameter of high amplitude points at time t 
(DIAMfiait) denotes the diameter of high amplitude points). 

We also keep track of where the walk was started. Since we want to be able to compare 
values on the two different domains, we will need to compute relative rather than actual 
values: Instead of listing a™"(l) and S{1), we list ^.e^l) = 7^^) and Srei{l) = ^^S(s)- 

The following six figures, numbered 15-20, show some typical plots obtained from our nu- 
merical calculation. We would like to point out that there is a lot of variation in the character 
of the plots: On both domains, the plots of the four quantities look different for diff'erent 
starting nodes but in almost all cases, the four quantities approach an asymptotic limit after 
about 200 time steps. 



Figure 15: Square Snowflake: Random Walk Simulation 
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Figure 16: Square Snowflake: Random Walk Simulation 
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Figure 17: Square Snowflake: Random Walk Frame 




Figure 18: Triadic Snowflake: Random Walk Simulation 




Figure Uni shows a simulations that takes more than 200 time steps for the values to approach 
the steady state. 



Let us look at some random walks started at low entropy/low nodes on both domains. 



Table 3: Square Snowflake decreases) 







PR{oo) 


DIAM{oo) 






DIAMrJioo) 


0.4572 


8.21 * 10-3 


25.6456 


0.9254 


0.0559 


0.707 


0.5829 


0.3216 


0.12 


25.5995 


0.9512 


0.0667 


0.69 


0.54 


0.3124 


0.146 


25.5997 


0.9071 


0.0557 


0.6972 


0.6225 


0.156 


0.4777 


25.5996 


1.3436 


0.0541 


0.6972 


0.3871 


8 * 10-^ 


0.4777 


25.6527 


1.3295 


0.0562 


0.707 


0.3716 


9.7 * IQ-'^ 


0.4777 


25.6455 


1.5202 


0.0493 


0.7071 


0.0783 



The diameter increases with decreasing ()^^'*(1) values but the relative diameter of high 
amplitude points decreases. Note that the entropy of high amplitude points is lying at 
around 70% of the maximum value which is much larger than we would expect, but the 
minimum ratio distances are very low (~ 5% of the maximum value) . 

Also note the large diameter but very low relative diameter of the high amplitude points in 
the last row. The corresponding plots for this simulation are given in figure 21. 



Figure 21: Square Snowflake: Random Walk Simulation node 2415 
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At a time step of about 200, and S]^^ start jumping between two values; an indication 
of localization? 



Table 4: Square Snowflake {S{1) decreases) 







PR(oo) 


DIAM{oo) 






DIAMj;f{oo) 


0.314 


0.1464 


25.5997 


0.9071 


0.0557 


0.6972 


0.6225 


0.3436 


0.07 


25.6456 


0.888 


0.0566 


0.7071 


0.6483 


0.3572 


0.058 


25.5997 


1.1761 


0.0575 


0.6972 


0.2941 


0.4572 


8.21 * 10-3 


25.6456 


0.9254 


0.0559 


0.707 


0.5829 


0.492 


1.78* 10-3 


25.6456 


0.9135 


0.0559 


0.7071 


0.6047 


0.492 


1.52 * 10-3 


25.6005 


0.9209 


0.0568 


0.6972 


0.5998 



The diameter of high amplitude points is large and does not decrease with decreasing entropy 
values. Here, row 3 sticks out with a relatively low DIAMJ^^ [oo) value. The corresponding 
plots are given in figure 22. 



Figure 22: Square Snowfiake: Random Walk Simulation node 824 
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In both tables, £)^^'(oo) and Sf^^{oo) show little fiuctuations. 



Table 5: Triadic Snowflake decreases) 







PR{oo) 


DIAM{oo) 






DIAM^Jioo) 


0.4776 


0.063 


156.8692 


0.5492 


0.0153 


0.4897 


0.7333 


0.3973 


0.231 


169.2816 


0.7011 


0.0144 


0.4883 


0.5033 


0.0206 


0.4579 


221.154 


0.9773 


0.0231 


0.5247 


0.1822 


0.0163 


0.4663 


86.536 


0.011 


0.0147 


0.4897 


0.1895 


0.0617 


0.4285 


221.3947 


0.9558 


0.023 


0.5246 


0.1944 


1.13* 10-=^ 


0.4201 


222.6706 


1.0908 


0.0308 


0.5301 


0.187 



In contrast with what was observed for the square domain, we see here in table 5 different 
participation ratios for different starting nodes. Note that on the triadic domain the entropy 
of high amplitude points are lying at only about 50%, as opposed to about 70% on the other 
domain. As observed on the square domain, DIAMJ^'^'- drops with decreasing c?™™(l). This 
is not the case when we decrease S{1), as shown in table 6. 



Table 6: Triadic Snowflake (5'(1) decreases) 







PR{oo) 


DIAM{oo) 


^ra(oo) 




DIAMll^ioo) 


6.25* 10-^ 


0.4277 


214.5546 


1.112 


0.0348 


0.5263 


0.1964 


0.3973 


0.231 


169.2816 


0.7011 


0.0144 


0.4883 


0.5033 


0.4776 


0.063 


156.8692 


0.5492 


0.0153 


0.4897 


0.7333 


0.4752 


0.063 


319.6411 


0.5881 


0.0163 


0.4897 


0.7852 


0.4895 


0.05 


222.6769 


0.5999 


0.025 


0.5310 


0.6707 


0.4969 


0.042 


222.6763 


0.6102 


0.025 


0.531 


0.6707 



Unlike in the case of the square Snowflake, we notice here that a decrease in the initial node 
entropy results in a decrease in DIAM{oo). 



We next present in tables 7 and 8 some random walks started at high entropy/high ^^j" 
nodes on both domains. We point out that on both domains, no node with i9™j"(l) and S{1) 
greater than 1/2 of their respective maximum values exists. 



Table 7: Square Snowflake 







PR{oo) 


DIAM{oo) 






DIAMJl^Jioo) 


0.7109 


0.3328 


25.6465 


1.2798 


0.0644 


0.707 


0.2032 


0.7108 


0.3312 


25.5997 


1.2847 


0.065 


0.6972 


0.1961 


0.5625 


0.3394 


25.762 


0.9719 


0.0548 


0.7071 


0.5256 


0.5624 


0.3375 


25.7173 


0.9733 


0.0563 


0.6972 


0.5315 


0.3332 


0.5414 


26.6456 


1.2564 


0.0538 


0.7071 


0.2661 


0.3125 


0.4227 


25.5996 


1.0088 


0.0796 


0.6972 


0.4651 



Entropy values of high amplitude points are the same as for starting nodes with low 9™j"(l) 
and S{1) values. DIAMJ^l'' {oo) varies a lot. 
Table 8: Triadic Snowflake 





SrelW 


PR{oo) 


DIAM{oo) 






DIAMl^^ioo) 


0.575 


0.0113 


211.4606 


0.8316 


0.0229 


0.5256 


0.275 


0.527 


0.3067 


222.0036 


0.7875 


0.0236 


0.5265 


0.3438 


0.4994 


0.4138 


222.0426 


0.8087 


0.0236 


0.5265 


0.3155 


0.4598 


0.5141 


222.0645 


0.8299 


0.0236 


0.5265 


0.2896 


0.3065 


0.743 


86.2551 


0.9098 


0.0151 


0.4939 


0.1511 


0.3248 


0.5546 


87.0524 


0.8781 


0.026 


0.4896 


0.1892 



DIAMJ^l^ (oo) values are lower than on the square domain. DI AM (oo) also seems to be 
lower and D)^^'(oo) only lies at about 2% of the maximum value instead of at about 5.5% 
on the other domain. Again, participation ratios vary on the triadic domain while they are 
constant for the square Snowflake. 



Tables 9 and 10 refer to the simulations run with high and low S{1) values. 



Table 9: Square Snowflake 







PR(oo) 


DIAM{oo) 






DIAMlfioo) 


0.7656 


1.31 * 10-^ 


25.5997 


1.2649 


0.0683 


0.6972 


0.2045 


0.9376 


6.49 * 10-5 


25.6468 


0.8764 


0.0548 


0.7071 


0.6714 


0.9376 


2.96 * 10-^ 


25.6473 


0.9209 


0.0592 


0.7071 


0.596 


0.9376 


2.1 * 10-9 


25.6053 


0.8643 


0.0501 


0.6972 


0.702 


0.2188 


4.07* 10-12 


25.6595 


0.925 


0.0533 


0.7071 


0.5922 


0.9376 


6.19 * 10-^ 


25.5997 


0.8796 


0.0563 


0.6972 


0.6725 



Table 10: Triadic Snowflake 







PR{oo) 


DIAM{oo) 






DIAMZ^ioo) 


0.4952 


0.04 


222.6771 


0.6844 


0.025 


0.531 


0.4638 


0.59 


4.32 * 10-=^ 


221.6813 


0.6529 


0.0182 


0.5228 


0.5147 


0.6383 


0.6554 


222.6791 


0.7485 


0.0253 


0.5324 


0.3408 


0.734 


0.147 


210.968 


0.747 


0.0219 


0.524 


0.3908 


0.772 


0.1176 


213.7942 


0.7185 


0.0224 


0.5255 


0.4135 


0.8483 


1.05* 10-6 


211.4231 


0.717 


0.0229 


0.5256 


0.3808 



Note that low S{1) values are not sufficient for DIAMf^f{oo) values to drop. 



Finally, tables 11 and 12 provide the results of the simulations run with low ^^"(1) and 
high S{1) values. 



Table 11: Square Snowflake 



<9re/(l) 




PR{oo) 


DIAM{oo) 






DIAM^Jioo) 


0.1876 


0.9554 


25.5997 


1.3179 


0.0624 


0.697 


0.1221 


0.1388 


0.949 


25.6047 


0.9904 


0.0563 


0.6972 


0.5015 


0.0936 


0.9363 


25.5997 


1.3681 


0.0692 


0.6972 


0.087 


0.09 


0.9358 


25.6455 


1.3514 


0.0548 


0.7071 


0.0842 


0.078 


0.9044 


25.6456 


1.3486 


0.0548 


0.7071 


0.0882 


0.078 


0.9044 


25.5996 


1.3208 


0.0579 


0.6972 


0.1165 



Note the high and low values of the relative diameter of high amplitude points of row 2 an 
3, respectively. 



Table 12: Triadic Snowflake 





SrelW 


PR{oo) 


DIAM{oo) 






DIAMl^^ioo) 


0.206 


0.7478 


86.5501 


0.8474 


0.0292 


0.4897 


0.2385 


0.1618 


0.9957 


86.6407 


0.8497 


0.0313 


0.4897 


0.2546 


0.1147 


0.8944 


222.1414 


0.8283 


0.0238 


0.5282 


0.3386 


0.09 


0.8361 


213.3961 


0.8561 


0.0236 


0.5254 


0.1841 


0.082 


0.7879 


86.8069 


0.8565 


0.0312 


0.4897 


0.275 


0.061 


0.7016 


213.8787 


0.9734 


0.0233 


0.5243 


0.121 



Even though we start in a high entropy region, it seems to be sufficient to have low drei{^) 
for a drop in DIAMJl^^ (oo) to occur. 



7 Results and Discussion 



Our numerical calculations show the following results on both domains: 

• All five quantities of interest approach an asymptotic limit. 

• D'^^^{oo) and Sj^^{oo) are independent of the starting node. Values lie at about 5.5% 
and 70% of the respective maximum values on the square Snowfiake and at about 2.5% 
and 50% of the respective maximum values on the triadic Snowfiake. 

• The asymptotic limit of the Diameter {DI AM{oo)) varies a lot within a region but 
overall seems to be a little lower when walks are started in low distance/low entropy 
regions. 

• When 5'(1) and are very low, the diameter of high amplitude points (D/AM^^'(oo)) 
is significantly smaller than for walks started in any other region. 

• No node with 5'(1) and d'^^^{l) values larger than half their respective maximum values 
can be found. 

Comparing results on the triadic and square domains, we observe the following: 

• Astonishingly, the participation ratio of the square Snowfiake is approximately the 
same for any starting node. On the triadic domain, we find three different approximate 
values: even if S{1) and c}™j"(l) values are similar (row three and four in table Ej), we 
can observe different participation ratios. 

• The participation ratio on the triadic Snowfiake is higher than on the square Snowfiake: 
On the square domain the participation ratio lies at 0.5%, the three values on the 
triadic domain correspond to about 1, 3 and 4% of the total number of nodes in the 
discretization. 

• On the triadic domain, no node with 9™^" < |max(9™j") and S < ^max(S') can be 
found, while on the square domain, we can find a node satisfying S < max (5) for 

^rat ^ 2 y^rat ) ■ 



• The minimum ratio distances of high entropy points on the square domain are constant 
and lie at about 5% of its maximum value on the domain. On the triadic Snowflake, 
minimum ratio distances of high entropy points vary: we find values from 1% to 3% of 
the maximum. These seem to be independent of where the starting node lies. 

• For both domains, the S'^^^(oo) values do not show sensitivity to initial conditions. It 
is rather surprising, though, that they are larger on the square Snowflake. 

In summary, we can say that random walks started at low 9^^" nodes show localization in 
the form of decreased diameter of high amplitude points. This localization occurs on both 
domains with the same frequency. The value of the initial node entropy does not seem to 
have a great influence on the dynamics of the system: a decrease does not necessarily result 
in a decrease of the diameter of high amplitude points. DJ^^J^(oo) and Sl^{oo) always level 
off to about the same value on both domains, with the (oo) very low and Sl'^{oo) values 
relatively large. This is a little surprising and requires further investigation. The main 
result of the numerical simulations is that the dynamics of the two systems are very similar. 
Hence, we can not conclude that narrow channels in the geometry of the square domain cause 
localization of the Dirichlet eigenfunctions. But there certainly are noticeable geometric 
differences on the two domains. The nodal entropies found on the square Snowflake are 
signiflcantly higher. Maybe the greater diversity in distribution on this domain is the 

cause for localization found in some of the Dirichlet eigenfunctions. This hypothesis should 
be further investigated. One should create weighted graphs with different nodal entropy 
distributions and test how the dynamics of the system is influenced. Another direction of 
further research in this area would be to conduct a thorough analysis of the transfer matrix 
T. It would be interesting to determine if the nature of the eigenvalues and eigenfunctions 
of T reveals information about the course taken by a random walker. 
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